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The boundary condition dependence of the critical behavior 
for the three dimensional Anderson transition is investigated. 
A strong dependence of the scaling function and the critical 
conductance distribution on the boundary conditions is found, 
while the critical disorder and critical exponent are found to 
be independent of the boundary conditions. 



The Anderson transition is a continuous quantum phase 
transition which separates metallic and insulating phases of 
the non- interacting electron gas [Q,^ . It is expected that the 
critical behavior should depend only on the basic symmetry 
of the Hamiltonian under the operation of time reversal . 
Recently this was clearly confirmed by numerical simulation 

i- 

At the transition the correlation length diverges and quan- 
tities such as level statistics ^||^ and the conductance distri- 
bution ]5p[-p^ become size independent and universal. The 
discovery that the critical level spacing distribution 
and indications that the critical conductance distribution 
[ p^lJl^ also, depend on the boundary conditions was unex- 
pected. 

In this letter, we analyze both the corrections to scaling in- 
duced by different boundary conditions and the effect of the 
different boundary conditions on the critical disorder, critical 
exponent and critical conductance distribution. None of the 
boundary conditions we consider break time reversal symme- 
try, so no change in the critical behavior can be predicted 
on the general grounds of a transition between universality 
classes. 

We report the results of two different simulations. In the 
first the scaling behavior of the localization length of electrons 
on long quasi-lD bars is examined. In the second the con- 
ductance distributions for ensembles of cubes of disordered 
material in a two probe measuring geometry are determined. 
Both simulations were repeated for three different bound- 
ary conditions; pbc) periodic boundary conditions in both 
transverse directions, mbc) periodic boundary conditions in 
one direction and fixed boundary conditions in the other and 
fbc) fixed boundary conditions in both transverse directions. 
These boundary conditions can be thought of as correspond- 
ing to different topologies; fbc corresponds to the topology of 
a wire, while mbc corresponds to that of a hollow cylinder. 

We find that the location of the mobility edge separating 



the localized and diffusive phases is unaffected by the choice 
of boundary condition. This is also true for the critical expo- 
nent. However, the scaling function of the localization length 
and the critical conductance distribution are found to depend 
strongly on the choice of boundary condition and hence on 
the topology of the sample. 

For the numerical simulations we have used the Anderson 
model 
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where C\(Ci) denotes the creation (annihilation) operator of 
an electron at the site i of a 3D cubic lattice. Energies Wi de- 
note the random potential distributed independently and uni- 
formly in the range [— V(^/2]. The hopping is restricted 
to nearest neighbors and its amplitude is assumed to be the 
energy unit, V = \. 

We consider first a quasi-ld system, i.e., a long bar, with 
cross section L x L. The standard transfer matrix tech- 
nique allows us to calculate the localization length of electrons 
X{Ef, W, L) on the bar within a desired accuracy p^ , p^ . The 
dependence of the quantity 
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on the width of the bar is then analyzed using the finite size 
scaling method. In this simulation we set the Fermi energy at 
the band center Ef = and vary the strength of the random 
potential W and the cross section size L. 

The data obtained are shown in Fig.|l| In the absence of 
any corrections to scaling, plotting A vs W should show the 
critical disorder Wc as the common crossing point of the data. 
However, as seen in Fig.^, the curves for different sizes do not 
cross at a common point. For pbc a previous analysis ^] 
has suggested that main reason for this is the existence of a 
correction due to an irrelevant scaling variable. For mbc and 
fbc surface and edge contributions might also be important. 
Such corrections are irrelevant. For example, for a surface 
effect we expect the corrections vanish as L~^. 

To take account of corrections to scaling we use the method 
described in [|l^. The data are fitted to the scaling form 



A = F(V)L'/",0L^) 
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where v is the critical exponent describing the divergence of 
the localization length, (p is the leading irrelevant variable and 
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y its irrelevant exponent. The best fit is found by minimizing 
the statistic in the usual way. 

When fitting the data, Eq.(^ is expanded in a Taylor series 
to first order in the irrelevant variable. 



(4) 



Both Fo and Fi are expanded to third order in their argu- 
ments. The relevant and irrelevant scaling variables are ex- 
panded in power series of the dimensionless disorder w — 
(Wc - W)/Wc as follows 



tp = tplW , (j> = 4>o 
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For pbc the expansion of the relevant field in (|5|) was continued 
to quadratic order as this gave a better quality of fit. Also, 
for fbc, data for L — 4 had to be omitted in order to obtain 
an acceptable fit. 

The numerical data and the associated fits are shown in 
Figs. |l|j^ and ^ The estimates for the critical parameters 
are listed in Table ^ and further details of the fits in Table ^ 
It can be seen from Fig. ^ that the main effect of imposing 
a fixed boundary condition at a given disorder is a decrease 
of the localization length A. Also the deviations from scaling 
behavior are much larger than those for periodic boundary 
conditions. For fbc the corrections are around 15% for the 
smallest system size, while they are only around 2% for pbc. 
The estimate of the irrelevant exponent for pbc is consistent 
with that in jl7|. The estimates for mbc and pbc are both 
much larger and close toy = —1 suggesting that the dominant 
correction is a surface effect. 

In Fig. ^ we see that taking account of corrections to scal- 
ing by plotting Acorrcctcd = A — cj>L^ Fi(4>L^^'^) as the ordinate 
restores a common crossing point to the data for different sys- 
tem sizes for each boundary condition. Further the critical 
disorder seems to be the same for all three boundary condi- 
tions. This is reinforced by looking at Table | where 95% 
confidence intervals are also given. Thus the location of the 
mobility edge does not seem to be affected by the choice of 
boundary conditions. The same is true for the critical expo- 
nent. 

However, the scaling function, and in particular Ac — 
Fo{0), do depend strongly on the boundary conditions. The 
estimates of Ac given in Table | are widely separated with 
no overlap of the confidence intervals. The change in the es- 
timated Ac is roughly 25%, while the analysis suggests that 
Ac has been estimated to within an accuracy of about 5% for 
fbc and to less than 1% for pbc. (Note that the apparent 
differences in the scaling function in the localized regime in 
Fig. 1^ are not important; these are an artifact of the scaling 
procedure in which the absolute value of the correlation and 
localization lengths are not determined.) 

Next we look at the critical conductance distribution Pc(ff). 
Since the correlation length diverges at the critical point, the 
critical conductance distribution of a phase coherent conduc- 
tor should be scale invariant [^. This was confirmed in nu- 
merical simulations |^,^ and it was also confirmed that Pc{g) 
depends on the universality class The conductance of a 
classical conductor depends on the aspect ratio (i.e. the ratio 
of the cross section to the length) but not on the shape of 
its cross section or on its topology (i.e. whether it is a bar 



or a cylinder.) Thus, a dependence of the critical conduc- 
tance distribution of a quantum conductor on its aspect ratio 
is also expected. Whether or not the the critical conductance 
distribution should depend on the conductor's cross sectional 
shape and its topology (boundary conditions) is less clear. 

We simulated the conductance distribution for an ensemble 
oi L X L X L cubic samples in a two probe measuring geom- 
etry using a Green's function iteration technique uM. The 
conductance in units of /h is 



g = 2tr(ttt) 
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where t is the transmission matrix found in the Green's func- 
tion iteration and the factor of two takes account of spin de- 
generacy. 

We first simulated an ensembles of 1,000,000 systems for 
each boundary condition. The resulting distributions are 
shown in Fig. ^. The choice of boundary condition affects 
Pc{g) especially for small g. There is a tendency towards 
more insulating behavior for fbc. 

It is expected on general grounds that Pc[g) will be size 
independent for "large enough" system sizes. This leaves open 
the possibility that the dependence of Pc{g) on the boundary 
conditions in Fig. ^ may simply be an indication that the 
system sizes employed in our work are not "large enough". 
We therefore decided to examine the size dependence of Pc(g). 
The dependence of < g > on the size L for L = 4 to L = 20 
for each boundary condition is shown in Fig. Ensembles of 
100, 000 systems were generated for each system size with the 
exceptions of L = 16 and L = 20 where the ensemble sizes 
were reduced to 25, 000 and 10, 000 respectively. 

To estimate the asymptotic value of < </ > in the limit 
L ^ cxD we have assumed that the size dependence is due to 
an irrelevant scaling variable and fitted the data to 



< g{L) > = < g{(x) > +aL^ 
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The details of the fits are given in Table III . From the results 
it seems clear that the boundary condition dependence will 
not disappear as L cxd. The asymptotic value of the mean 
conductance is reduced by about 40% while the analysis sug- 
gests that it has been estimated to within a few percent. The 
numerical data are consistent with a size independent and 
universal Pc{g) but also one that depends on the boundary 
conditions even in the limit L oo. 

The size dependence of the critical conductance fiuctua- 
tions were analyzed in a similar way and the results are also 



given in Table [II, (For fbc no size dependence of the vari- 
ance was detected within the accuracy of the simulation so 
the value given is simply a weighted average over different 
system sizes.) The values given can be compared with the 
value of var((7) = 1.18 (with g in units of e^/h) for universal 
conductance fiuctuations in the metallic regime [^. Though 
the absolute magnitude of the fiuctuations is smaller than in 
the metallic regime, the conductance fluctuations are in fact 
of the same order of magnitude as < g > at the transition. 
The estimates of the irrelevant exponent in Table 
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again 

suggest the presence of surface effects. However, for the con- 
ductance calculation it seems to be pbc which presents a sur- 
face effect rather than fbc. For the latter boundary condition 
a surface correction does not seem to be present and we find a 
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value of the irrelevant exponent which is consistent with that 
for the scaling behavior of the localization length with pbc 
given in Table |. 

In conclusion we have found that the some aspects of the 
critical behavior at the Anderson transition, the location of 
the mobility edge and the critical exponent, are independent 
of the boundary conditions. Other aspects, the scaling func- 
tion and the critical conductance distribution, seem to depend 
strongly on the boundary conditions. The invariance of the 
critical disorder is reasonable since localized states should be 
unaffected by boundary conditions. However, this argument 
does not hold in the critical regime and, as we have shown, 
important aspects of the critical behavior do depend on the 
topology. The exception is the critical exponent provided, as 
here, that the boundary conditions do not break time reversal 
symmetry. 

Finally, we recall that a similar dependence of the scal- 
ing functions on boundary conditions has also been found in 
classical percolation [ po| , pl| . This suggests that the quantum 
nature of the Anderson transition might not be crucial for the 
existence of this sort of effect. 

The authors would like to thank K.Fukushima, S.Todo, 
Y.Okabe and N.Hatano for useful discussions. 
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pbc 
mbc 
fbc 


1.56(55,58) 
1.60(56,64) 
1.54(41,61) 


16.54(53,55) 
16.47(42,52) 
16.49(39,64) 


0.576(.574,.577) 
0.502(.494,.509) 
0.426(.403,.442) 


-2.8(3.2,2.4) 
-1.3(1.5,1.2) 
-1.2(1.4,1.0) 



TABLE I. The best fit estimates of the critical exponent, 
the critical disorder, Ac and the irrelevant exponent together 
with their 95% confidence intervals. 











Q 


pbc 


336 


12 


335 


0.3 


mbc 


238 


11 


233 


0.4 


fbc 


204 


11 


212 


0.2 



TABLE IL The boundary condition, the number of data 
the number of parameters A'p, the value of for the 
best fit and goodness of fit Q. The system sizes used were 
1/ = 4, 5, 6, 8, 10, 12, 14 with the exception of fbc where the 
data for L = 4 were omitted and pbc where data for L = 16 
are also included. The accuracy of the numerical data was 
either 0.1% or 0.05%. The range of disorder used was W = Vb 
to W = 18. 
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< ff(oo) > Q y' 


var(g(cxj)) 


pbc 0.89 ±.02 0.2 -0.9 ±.1 
mbc 0.71 ±.01 0.6 -1.1 ± .1 
fbc 0.560 ±.002 0.9 -2.3 ± .7 


0.472 ± .01 
0.408 ± .003 
0.349 ± .001 



TABLE III. The estimated value of < 3 > in the ther- 
modynamic limit L — » 00 together with standard errors, the 
goodness of fit Q and an estimate of the irrelevant exponent. 
An estimate of Ya.r{g) in the limit L ^ 00 together with stan- 
dard errors is given in the last column. 
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FIG. 1. A vs. W. for periodic (upper curves), mixed (mid- 
dle curves) and fixed (lower curves) boundary conditions. 



FIG. 4. The dependence of the critical conductance dis- 
tribution on the choice of boundary conditions. Here the 
Fermi energy Ef = 0,the system size L = 10 and the disorder 
W = 16.54 independent of the choice of boundary conditions. 



(1)0.6 



< 0.5 




0.7 



Si . 6 




FIG. 2. A vs. W after the surface corrections are removed. 
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FIG. 5. The size dependence oi < g > for different bound- 
ary conditions. Here the Fermi energy Ep = 0.5 for which we 
estimated Wc = 16.53 independent of the choice of boundary 
conditions. 
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FIG. 3. The s caling functions for -different boundary con- 
ditions. 



4 



